Experimental confirmation of tissue liquidity based on the exact solution of the Laplace equation 
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The notion of tissue surface tension has provided a physical understanding of morphogenetic phenomena 
such as tissue spreading or cell sorting. The measurement of tissue surface tension so far relied on strong ap- 
proximations on the geometric profile of a spherical droplet compressed between parallel plates. We solved the 
Laplace equation for this geometry and tested its solution on true liquids and embryonic tissue fragments as 
well as multicellular aggregates. The analytic solution provides the surface tension in terms of easily and ac- 
curately measurable geometric parameters. Experimental results show that the various tissues and multicellular 
aggregates studied here are incompressible and, similarly to true liquids, possess effective surface tensions that 
are independent of the magnitude of the compressive force and the volume of the droplet. 

PACS numbers: 68.03.Cd, 68.05.Gh, 87.17.-d, 87.18.La 



In the absence of external forces, an embryonic tissue frag- 
ment (typically less than 1 mm in size) rounds up similarly to 
a liquid drop, to minimize its surface energy. Cells of two dis- 
tinct tissues randomly intermixed within a single multicellu- 
lar aggregate sort into separate regions similarly to coalescing 
immiscible liquids. To account for these observations Stein- 
berg formulated the Differential Adhesion Hypothesis (DAH), 
which states that embryonic tissues or, more generally, tissues 
composed of motile cells with different homotypic adhesive 
strengths should behave analogously to immiscible liquids 
il, 2]. DAH implies that such tissues, similarly to ordinary 
liquids, possess measurable surface and interfacial tensions 
generated by adhesive and cohesive interactions between the 
component subunits (molecules in one case, cells in the other). 
Predictions of DAH have been confirmed both in vitro [3] and 
in vivo [3, El 0], and the surface tensions of different em- 
bryonic tissues were measured and the values accounted for 
the observed mutual sorting behavior (st [j^ . Currently, the 
only available method to measure the surface or interfacial 
tension G of submillimeter size droplets of tissue aggregates 
is by compression plate tensiometry [7, 8]. The method, based 
on the Laplace equation, so far relied on various approxima- 
tions of the geometrical profile of an equilibrated spherical 
droplet compressed between two parallel plates and yielded 
G values that can at best only be considered relative and their 
independence on droplet size and compressive force question- 
able. In this Letter, by analytically solving the corresponding 
Laplace equation, we determine the exact profile of a com- 
pressed droplet, which allows, for the first time, to accurately 
and reliably determine the absolute value of tissue surface ten- 
sions in terms of easily and accurately measurable geometric 
parameters. Furthermore, we show that the method can read- 
ily be extended to the simultaneous compression of several 
droplets of different sizes. Finally, we apply our new method 
to a large number of compression plate measurements for cal- 
culating G of true liquids (for validation) and several tissue 
and multicellular aggregates. Our results show that the stud- 
ied systems are incompressible and the obtained surface ten- 



sions are independent of the magnitude of the compressive 
force and the volume of the droplets. Our results provide 
strong evidence for the concept of embryonic tissue liquidity 
and support its usefulness for the interpretation of early mor- 
phogenetic processes. More importantly, since surface tension 
is a measure of the liquid's cohesivity, in the case of tissues it 
must be related to molecular parameters. Indeed, it was shown 
on theoretical grounds that G ^ JNz, where /, N and T are re- 
spectively the bond energy between two homotypic cell adhe- 
sion molecules (CAMs), the surface density of CAMs and the 
effective life time of the adhesive bond [8]. The linear depen- 
dence of CT on has recently been confirmed experimentally 
. Thus our method of determining the absolute value of G 
has important biological implications as it quantitatively re- 
lates a macroscopic tissue property to biomolecular entities. 

The surface tension O" of a small liquid droplet compressed 
between two pressure plates can be determined from its geo- 
metric shape (Fig.[T]\-B). Here we consider droplets with ra- 
dius Rq much smaller than the corresponding capillary length 
Rc ^ (cr/pg)^/^, for which the effect of gravity can be ne- 
glected (e.g., for water, with G = 0.07 N/m and density p = 
lO^kg/m^, Rc ~ 2.7 mm). Thus, the shape of a submillime- 
ter liquid drop placed on a horizontal plate (Fig. [TJ\) is a 
spherical cap of radius Rio and height Hq. With these pa- 
rameters, simple geometric considerations provide: (i) the 
(complementary) contact angle 6 = cos~^ {Hq/Riq — 1), and, 
assuming incompressibility, (ii) the radius of the suspended 
drop, Rq = /?io[(2 - cose)cos'^(e/2)]. While Hq and Rio 
can be measured with high accuracy (e.g., < 1%), the rela- 
tive error AG/G ^ [(1 + cos G)/G sin G] {AHq/Hq + ARiq/Riq) 
can still be very large (e.g., - 30% for G = 20° and > 130% 
for G < 10°). Thus, for determining the surface tension it is 
desirable to reduce the adherence between the drop and plates 
and avoid using quantities that explicitly contain the contact 
angle. 

The compressed drop (Fig. [T^) has rotational symmetry 
about the z-axis and reflection symmetry with respect to the 
equatorial plane z = H/2.ln this plane the surface of the drop 
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FIG. 1: Diagram of an (A) uncompressed and (B) compressed liquid 
drop. (C) Snapshots of uncompressed (top) and compressed (bottom) 
cushion tissue (CT) droplet in culture medium. (D) Same for two 
different-size water drops in olive oil. 



pressed liquid droplet (Fig.[T^). This can be obtained by the 
exact integration of the Laplace equation 



Ap = G 



_r(l+z'2)i/2 (l+z'2)V2 

subject to the boundary conditions (see Fig.[T^) 

z{R3)=0, z'(7?3)=tane, 
z{Ri)=H/2, z\Ri)=oo. 



— const, 



(2) 



(3) 



The terms in the square brackets in Eq. © represent the prin- 
cipal curvatures of the drop's surface at a point determined by 
zir). 

We are now in the position to determine cr in terms of 
the easily and accurately measurable quantities H, Ri and F. 
Equations ([B imply that both R2 and /?3 can be expressed in 
terms of Ri 

R2=Ri/{2a-l), R3=pRi, (4a) 
where the dimensionless parameters a and j3 are given by 

a = Ap/{2a/R,), (4b) 



has two principal radii of curvatures Ri and R2 (Fig.[T^). R3 
is the radius of the droplet' s circular area of contact with ei- 
ther compression plates. The degree of compression depends 
on the magnitude of the compression force F applied to the 
upper (or lower) plate. In terms of Ri and R2 the excess pres- 
sure inside the drop due to the surface tension is given by the 
Laplace formula Ap = g{1 / Ri -\- 1 / R2) . Thus, at mechanical 
equilibrium, the balance of forces dictates the following rela- 
tions, respectively valid at the upper (or lower) plate and the 
equator 



F = Ap tzR] - 27r7?3a sin 

= 7tG[Rl{\/Ri^\/R2)-2R3^me] , 

F = Ap 7iR\ -2nRia = noRi {R1/R2-I) 



(la) 



(lb) 



Although H, Ri and F can be easily measured, the accurate 
determination of a directly from either of these two equations 
also requires the quantities R2, and R3 that can only be mea- 
sured with large errors (especially the latter two). This prob- 
lem is often circumvented by assuming = (no adhesion to 
the plates) and/or making approximations on the lateral profile 
of the drop, e.g., R3 = Ri — R2 (the profile is a semicircle) or 
Ra=Ri-R2^[Rl- {H/2f] (the profile is a circular arc) 
[Sl ISO. Each of these schemes fails to give consistent results 
in some range of the compressive force or contact angle. For 
example approximating the lateral profile with a circular arc 
implies that the contact angle depends on the magnitude of the 
compressive force, a physically absurd conclusion. 

Clearly, the precise determination of the radii R2 and R3 
requires the knowledge of the exact profile z{r) of the com- 



i3=i3e(a) = (2a)- 



sine + A/sin^e+4a(a- 1) 



. (4c) 



Integrating Eq. Q with the boundary conditions ^ leads to 
an implicit equation for a 



i 

^=/e(«) = j z'{x)dx 



z!{x) 



1 



a 



1 



-1/2 



(5) 



and to the lateral profile of the compressed drop 

rr/Ri 



z{r) =Ri r \\x)dx. 
Jl5 



(6) 



Because the value of the integral in Eq. (0 depends weakly on 
the lower integration limit j3, one finds that a, and therefore 
7?2, are relatively insensitive to even large (~ 10°) sample- 
to-sample fluctuations of 6 (e.g., caused by local inhomo- 
geneities or impurities). By contrast, since j3 itself strongly 
depends on 0, so does R3. This explains why its determina- 
tion by fitting the lateral profile of the drop by circular arcs is 
impractical. 

Finally, cr can be expressed from either Eqs. ([T]), which in 
terms of a ^ f^^ {H/2R\) [see Eq. ([5])] can be rewritten as 

Ao = i^/27r7?iCJo = (a/cJo)(a-l) , (To = lmN/m. (7) 

This result provides a simple recipe to evaluate a from the 
measurement of H,R\ and F. The efficiency of the proposed 
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FIG. 2: Implementation of Eq. to evaluate the experimental data 
obtained in the compression measurements. The surface and inter- 
facial tensions are obtained from the slopes of the linear fit to the 
data points including the origin: (a) pure liquid (W-0, W-T), and (b) 
tissues and multicellular aggregates (CT, CHO, HUVEC, HUSMC). 



method can be further enhanced by simultaneously compress- 
ing several drops as shown in Fig.[TJ). For such compressions 
the quantities Aq and a in Eq. (|7]) need to be replaced respec- 
tively by F /IncoRi and Y^n^inOCn^ with = Y.n^\n and 

an=f^\H/2Rin)- 

It should be emphasized that Eq. ([5]) is valid for incom- 
pressible and compressible liquids alike. In addition, for in- 
compressible liquids volume conservation yields 



Jl5 

For 6 < 20°, to a very good approximation, ge{cc) 
Eliminating a between Eqs. (0) and ^ leads to 

H/2Ri = Ue{H/2Ro) ^ Uq{H/2Rq) , 



(8) 



(9) 



where is a universal function determined by the functions 
/o and go 1110(1 . Thus, when the adhesion between the drop 
and plates is weak (i.e., Q < 20°) there is a universal relation- 
ship between H/2R\ and H / 2Ro, valid for any incompressible 
liquid drop regardless of its type or size. 

We have tested the above theory to determine the interfa- 
cial tension of true immiscible liquids (water in olive oil (W- 
O) and turpentine (W-T)), where results obtained by other 
methods are available. Subsequently, we applied the ana- 
lytical results to determine the absolute values of tissue sur- 
face tensions. Measurements were performed using a specifi- 
cally designed compression tensiometer (for details see [11]). 
In the case of ordinary liquids, one or two spherical water 
drops of various size (ranging from 0.5 to 1 mm in diame- 
ter) were compressed simultaneously (Fig. [TJ)). In each ex- 
periment, drops were exposed to at least 3 successive com- 
pressions (up to 9) of increasing force and their profile was 
recorded at shape equilibrium. 



TABLE I: Surface and interfacial tensions of the studied systems ob- 
tained from the data shown in Fig. [3 The absolute errors Ac are 
standard deviations. The last two columns contain respectively the 
percentage relative errors and the number of data points. 
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FIG. 3: Implementation of the universal function given in Eq. ^ 
(solid curve) to assess the incompressibility of the (a) liquid and (b) 
tissue and multicellular droplets used in the compression measure- 
ments. 



For biological measurements we either used intact tissue 
(fragments of excised embryonic chicken cardiac cushions 
(CT) that round into spheres in about 12 hours) or model 
tissues (spherical aggregates composed of various cell types, 
such as Human Umbilical Smooth Muscle Cells (HUSMC), 
Human Umbilical Vein Endothelial Cells (HUVEC) and Chi- 
nese Hamster Ovary cells (CHO)). Multicellular aggregates 
were prepared as previously described flT]. Tissue droplets 
were compressed in culture medium at 37°C (Fig. [TJC). One 
to six droplets of each type (ranging from 250 to 600 /im in 
diameter) were compressed simultaneously. In order to avoid 
irreversible damage to the cells, no more than 2 compressions 
were performed on each droplet. 

According to the solution of the Laplace equation (|7]), for 
liquids, the pairs of data points {(a — 1), Aq} should lie on a 
straight line passing through the origin, with a given by the 
slope of the line. The data points are shown in Fig. [2l and the 
obtained results are summarized in Table [D Besides the values 
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of (7, the table also contains the standard deviation (Ad) of the 
measurements and the corresponding percentage relative er- 
rors (Act / a). The last column of the table lists for each system 
the total number of data points used and shown in Fig.O The 
obtained interfacial tensions of W-0 (18.6 ± 2.4 mN/m) and 
W-T (12.4 ± 1.5 mN/m) compare rather well with published 
data [12], i.e., 15.4 and 13.7 mN/m, respectively. It should be 
noted that the corresponding relatively large errors ( 12%) in 
Table IB reflect the sample sensitivity of G and not a deficiency 
of the method to determine it. Indeed, the large number of 
W-0 and W-T compression measurements (see Table |l|) were 
done in the course of four different days, each time using dif- 
ferent samples and compression plates. Also, most of these 
experiments were carried out using single droplets, which 
might have experienced somewhat different environmental 
conditions (sample-to-sample fluctuations). The actual errors 
Act corresponding to drops from the same batch were much 
smaller (~ 1%). Another possible source of error is related to 
the extent of compressions. In general, results from weaker 
compressions have larger errors. Indeed, the relative error of 
the surface tension Aa/a = AF /F ARi/Ri + Aa/(a - 1) 
[obtained from Eq. (|7])] may become very large for small com- 
pressions when F ^0 and a ^ I. 

By a similar procedure, the surface tensions of tissue 
droplets (i.e., interfacial tension in tissue culture medium) 
were also determined with the results shown in Fig. [2]3 and 
Table m The fact that the data points in Fig.[2]3 corresponding 
to each type of cell aggregates fall on well defined lines that 
go through the origin clearly demonstrates that embryonic tis- 
sues and cell aggregates have well defined surface tensions 
that are independent of their size or extent of compression. 
The small Act in case of CT (HUVEC) aggregates is due to 
the fact that the measurements were carried out using a single 
batch (of simultaneously prepared) aggregates. These aggre- 
gates were exposed to multiple multi-aggregate compressions: 
five or six aggregates compressed twice. In the case of CHO 
aggregates originating from a single batch, where a bigger 
number of compressions were performed on both single and 
multiple (two to five) aggregates, the relative error Aa/a was 
similar to that of W-0 and W-T. In the case of HUSMC four 
different batches of aggregates were used and only single ag- 
gregate compressions were performed. This may explain the 
relatively large error (~ 20%). Thus one may conclude that 
multi-aggregate compressions with droplets originating from 
a single batch is the most efficient and accurate way of exper- 
imentally determining tissue surface tension. This is partic- 
ularly important for highly viscous aggregates of living cells 
whose properties may change during their relatively long re- 
laxation time in a compression measurement. 

Furthermore, the fact that the data points {H /IRq^H /2R\} 



for all our measurements fall on the universal curve predicted 
by Eq. Q (Fig. [5]) demonstrates that the embryonic tissues 
and multicellular aggregates studied here, similarly to water, 
are incompressible. 

The results in Figs. [21 and [5] provide a wealth of infor- 
mation. They convincingly demonstrate that embryonic tis- 
sues or multicellular aggregates composed of embryonic cells 
(model tissues) indeed manifest liquid-like properties. Thus, 
such tissues are incompressible and can quantitatively be char- 
acterized in terms of liquid-like surface tensions. The mea- 
sured values of these parameters are independent of the size 
of the drops used or the magnitude of the force exerted on 
them. The employed method therefore provides the absolute 
values of tissue surface tensions. As discussed earlier, these 
results can be used to assess biomolecular quantities. 

The analytic solution of the Laplace equation presented 
here, combined with compression plate tensiometry, provides 
a novel, reliable and accurate way to determine liquid surface 
tension (and to our knowledge the only method applicable to 
tissues). Our approach may be particularly useful in the case 
when the interfacial tension of expensive liquids needs to be 
determined, where only small quantities of the materials are 
available. 

This work was supported by the National Science Founda- 
tion [FIBR-0526854]. 
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